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We obtain in a direct and rigorous manner a transition from a stable molecular hydrogen nH 2 
single chain to the quasiatomic two-chain 2nH state. We devise an original method composed of an 
exact diagonalization in the Fock space combined with an ab initio adjustment of the single-particle 
wave function in the correlated state. In this approach the well-known problem of double-counting 
the interparticle interaction does not arise at all. The transition is strongly discontinuous, and 
appears even for relatively short chains possible to tackle, n = 3-^6. The signature of the transition 
as a function of applied force is a discontinuous change of the equilibrium intramolecular distance. 

The corresponding change of the Hubbard ratio U/W reflects the Mott-Hubbard-transition aspect of 
the atomization. Universal feature of the transition relation to the Mott criterion for the insulator- 
metal transition is also noted. The role of the electron correlations is thus shown to be of fundamental 
significance in this case. The long-range nature of Coulomb interactions is included. 

PACS numbers: 31.15.A-, 71.27.+a, 67.80.F-, 67.63.Gh 


1. Motivation The metallization of solid hydrogen is 
one of the central problems in physics m, as well as 
in astrophysics of Jupiter, Saturn, and exoplanets [S]. 
The building block, the H 2 molecule, is relatively sim¬ 
ple, since the electrons are in the spin-singlet state in 
the ground state composed mainly of Is states of atoms. 
Those molecular states with the experimental value of 
bond length R « l.dOloo form at ambient pressure a 
molecular crystal with the lattice constant a ~ Tap S> R 
[7], where oq is the Bohr radius. Thus, to achieve atomic 
structure and metallicity one has to break the molecu¬ 
lar bond, e.g., by achieving a more typical atomic solid 
configuration with a R, what amounts in practice to 
applying an enormous pressure, as demonstrated repeat¬ 
edly over the decades mu- But even then one may not 
achieve metallicity, as due to the relatively large atomic 
separation the system may compose a Mott (or Mott- 
Hubbard) insulator |^I10|. as the original interelectronic 
correlation effects may still be sufficiently strong. The 
fundamental aspect of this paper is to address this issue 
in a rigorous manner, as detailed below. 

The related and probably more intriguing question re¬ 
lated to the hydrogen metallization is the conjecture that 
the system may exhibit a room-temperature supercon¬ 
ductivity. The line of reasoning HU EH is based on the 
Bardeen-Cooper-Schrieffer (BCS) theory that as the crit¬ 
ical temperature Ts is proportional to where M 

is the ionic mass, then the hydrogen metal should have 
Ts at least one order of magnitude higher than that of 
a typical metal. However, here again the correlation ef- 
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fects must be carefully taken into account as they should 
not be too strong to destroy the metallicity, but essen¬ 
tial, as they can lead by themselves also to the high- 
temperature superconductivity when the starting point 
is the atomic solid with half-filled Mott-Hubbard metal 
m- A marriage of strong electron-lattice and correla¬ 
tion effects should be accounted for at the starting stage. 
In this respect, here we determine rigorously the mag¬ 
nitude of the local correlation effects as described by 
the effective extended Hubbard model at the molecu¬ 
lar quasiatomic solid transition. The values of local 
electron-proton coupling have been calculated elsewhere 

US]- 

Recently, the molecular-atomic hydrogen transition 
has been discussed by a number of methods [TBUIH - In 
this respect, the present results can be regarded as a basis 
for testing the approximate methods. 

2. Model and method We start with the extended Hub¬ 
bard model with additional term Uion-ion expressing ion- 
ion repulsion namely. 
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where is the single-particle energy, tij are the so-called 
hopping integrals (all of them: to (intramolecular) and 
ti, and ^2 (intermolecular) have been marked explicitly in 
Fig.f^, U is the on-site Coulomb repulsion, and Kij is the 
amplitude of intersite Coulomb repulsion. The configura¬ 
tion of the molecules and related microscopic parameters 
are shown in Fig. Explicitly, those parameters can be 
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FIG. 1. Schematic representation of H 2 molecular chain with 
possible nearest-neighbor and next nearest neighbor hoppings 
{ti}i=o,i .2 marked. The labeling of the sites i = 1, 2, 3,... is 
specified, as well as the bond length R at intermolecular dis¬ 
tance a. The radius of included intersite Coulomb interaction 
is equal to 250a in the atomic representation, (see main text). 


defined in the following manner 

(wi| Hi \wj) = j (fr w*{Y)Hi{r)wj{v) 

= -\- tij{\ ) ( 2 ) 

{wiWj\Vi2\wiWj) = J d^rd^r' \wi{r)\‘^Vi2(r - r')\wj{r')\‘^ 

= C/^y+Ff„(l-%), (3) 

where Hi is the Hamiltonian of single electron in the 
medium and V 12 is the Coulomb repulsive interaction for 
a single pair of them. Furthermore, the Wannier func¬ 
tions are defined in terms of atomic (Slater) functions, 
i.e., 
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where the Is-type Slater function is defined as V’i(r) = 
ipir — Ri) = exp(—a|r — Ri|), with a being 

its inverse size in the medium, here taken as a variational 
parameter for given intermolecular distance a and deter¬ 
mined in the correlated state (i = R^ is the i-th proton 
position). Each Slater function is approximated by its 
expansion in the Gaussian basis 
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with p being the number of Gaussian functions taken 
into account to express accurately the functions {t/jiir)} 
and {Bg, Fg} is the set of adjustable parameters obtained 
from a separate procedure (for details see [H]). In effect, 
all the parameters: e^, Uj, U, Kij can be expressed in 
terms of a and are determined in the correlated state 
(after the exact diagonalization in the Fock space has 
also been carried out simultaneously). Strictly speaking, 


to obtain proper atomic limit we have to take: —)■ 

= G + (l/2)I]j(2/^*j + Kij) PO]. 

A brief methodological note is in place here. We di¬ 
agonalize exactly the clusters up to n = 5. However, 
to account for the long-range nature of Goulomb inter¬ 
actions, we include them in e^, tij, Kij, and the ion-ion 
interaction up to the distance 250a. In other words, we 
regard our cluser of nH 2 molecules as a single block of 
the length na immersed in the periodic system. In this 
manner, we utilize the periodic boundary conditions in 
the same sense as in the cluster perturbation approaches, 
cf. e.g. Ref. pT] . 

First, we determine values of hopping integrals to, ti, 
and ^2 Hubbard (intraatomic) and intersite interactions 
U and Kij, by determining the Wannier functions defin¬ 
ing them. The Wannier functions are expressed in terms 
of (atomic) Slater functions of adjustable size a~^, each 
of which is expressed in turn via p = 9 Gaussian func¬ 
tions m- In the next step, we determine the lowest 
eigenvalue of the Hamiltonian ([^ with periodic bound¬ 
ary conditions, by diagonalizing it by means of Lanczos 
algorithm. By finding the ground state energy in the 
Fock space and subsequently, by optimizing the energy 
also with respect to the size a~^, we obtain a true phys¬ 
ical ground state energy in the correlated state for given 
intermolecular distance a. Such extended calculations 
make the results not limited to the parametrized-model 
considerations, as U, tij, and Kij are evaluated explic¬ 
itly, as is also the ground state energy, for fixed a. One 
should note that in the procedure the bond length is also 
optimized (i? —>■ Reff)- All this is carried out first for the 
zero-applied force to obtain the system equilibrium con¬ 
figuration, i.e., the energy and the effective bond length 
Reff in the multimolecular configuration. In the second 
step, we apply the external compressing force / to the 
end molecules and determine the enthalpy minimum to 
trace the system evolution as a function of a. The part 
of the procedure connected with the optimization of the 
single-particle wave functions has been discussed in de¬ 
tail elsewhere m- The code used for the computation is 
also available |22| . 

Typical numerical procedure starts with fixing a and 
R, and selecting input value for the variational parame¬ 
ter a. Next, we vary the inverse wave function size a to 
find the physical ground-state energy employing, at every 
step, the so-called process-pool solution for the computa¬ 
tionally expensive problem of obtaining the microscopic 
parameters. We include all the intersite interactions for 
the sites with a relative distance smaller than 250a. The 
accuracy of the numerical results is set to be of the order 
of 10~'^Ry. Impact of these parameters was examined 
carefully in |20| and achieving the values above proved 
to be sufficient. The long range of the interactions in the 
atomic representation is expected to emulate correctly a 
longer-chain behavior. 

3. Results In Fig. we display the system enthalpy as 
a function of force / exerted on the molecules. We se¬ 
lected such anisotropic form of f, since a proper account 
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FIG. 2. Phase diagram encompassing quasiatomic and molec¬ 
ular states for nH 2 chain as a function of exerted force along 
the chain. Inset: Effective bond length Reff vs. the inter- 
molecular distance a. 


of hydrostatic-pressure case would require a considera¬ 
tion of square or 3-dimensional clusters, so the all direc¬ 
tions are physically equivalent. We assume a spatially 
homogeneous state of the system, in accordance with pe¬ 
riodic boundary conditions. Typically, two solutions ap¬ 
pear which we call respectively as the molecular and the 
quasi-atomic, each of which is characterized below. The 
solid points mark the transition points for 4i/2 (left) and 
3i?2 (right). In the Inset we show the effective bond 
length Reff vs a. Note two specific features. First, 
the lattice contracts in a discontinuous manner at the 
transition and second, the bond length jumps then from 
an almost single-molecule value Reff ^ a to the limit 
Reff ~ « signaling their separation into a quasi-atomic 
configuration. The last feature of the transition is ex¬ 
plicitly illustrated in Fig. where the projected electron 
density onto the plane composed of molecules is shown 
both on the molecular (top) and the quasi-atomic (bot¬ 
tom) sides of transition. The parameters a and Reff are 
listed also in each case. 

In Fig. 1^ the interdependence of the intermolecular dis¬ 
tance versus force is provided. As shown also in the 
inset to Fig. this figure demonstrates the first-order 
transition, since the cell volume a changes discontinu- 
ously at this mol. ^ quasiat. transition from 1.869ao to 
1.647ao. In the inset to Fig. we provide the ratio of 
the Hubbard interaction U to the bandwidth defined as 
W = —4fi-|-|fo + 2^2! —1^0 — 2f2|, as a function of /. Note 
the fundamental characteristic: this ratio jumps from the 
value U/W ~ 1.49 to 0.8. So, the system is indeed 
strongly correlated at the transition with the effective 
gap near the transition from the molecular side estimated 
as ig/W = U/W — 1 « 0.49. Furthermore, the value 
U/W ~ 1 is the typical value for the Mott-Hubbard tran¬ 
sition, weakly dependent on the system structure |23]. It 
is also interesting to note that here the Mott-Hubbard 
transition takes place from a non-magnetic (spin-singlet) 
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FIG. 3. Electron density projected onto the plane of H 2 
molecules near the molecular —> atomic transition: a) for 
the molecular state; b) for the quasiatomic state. In case 
a) the density is very small in the space in between molecules 
(vertical line joining the atoms). The critical intermolecular 
distance a and the effective bond length Reff are also speci¬ 
fied in each of the states. The corresponding critical force is 
fc = 25.705nA for 3772 system. 



FIG. 4. Equilibrium intermolecular distance a versus exerted 
force / for the case of 3772 system. The transition is of the 
first order as marked. Inset: U/W ratio vs. a showing that 
in the quasiatomic state U/W ~ 0.8 signaling the onset of a 
correlated metallic state. 


molecular configuration to a correlated atomic solid, with 
strong suggestion for metallicity {U/W < 1) of such sys¬ 
tem with one electron per atom. 

In Fig. 1^ we present a dependence of the intramolecu¬ 
lar (to) and intermolecular (ti) hopping parameters. At 
the transition the parameters to and ti roughly equal¬ 
ize showing again that the solid at small a is of quasi¬ 
atomic nature. In the inset we present the intramolecu- 
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TABLE I. Microscopic parameters at the molecular —> quasiatomic transition. The values are in eV if not specified otherwise. 



a (ao) 

R (ao) 

W 

to 


t2 

U 

Ko 

Ki 

K 2 

«(«o 

molecular 

1.869 

1.164 

17.881 

-15.487 

-6.161 

1.691 

26.472 

15.238 

13.516 

10.551 

1.194 

quasiatomic 

1.647 

2.386 

33.124 

-5.270 

-8.445 

0.164 

26.404 

10.825 

14.701 

8.976 

1.251 



FIG. 5. Evolution of the intra- and inter-molecular hop¬ 
ping parameters, to and t\ respectively, versus intermolecular 
distance a. Note the jump of a at the transition molecu¬ 
lar quasiatomic conhgurations as marked by the arrows. 
Inset: hopping probabilities: intra- and inter-molecular, Co 
and Cl, respectively. For discussion see main text. 


lar (Co = and intermolecular (Ci = 

hopping probabilities. In the molecular state Co ~ I and 
Cl « 0 (a strong molecular bond is formed), whereas 
in the quasiatomic state Cq Ci. These two quan¬ 
tities illustrate thus in a direct manner the nature of 
the states and in particular, the bonding in the quasi¬ 
atomic solid state. For the sake of completeness, in Ta¬ 
ble H] we list all the principal microscopic parameters in 
the correlated state for intermolecular distances at tran¬ 
sition. The values of Kij converge towards the asymp¬ 
totic value Kij —>■ 0 slowly. Such behavior is in accor¬ 
dance with the view [24] that the long-range part of the 
Coulomb interaction is not screened efficiently in low¬ 
dimensional systems. This is the reason why we have 
taken their long-range character. Parenthetically, note 
that we have neglected the direct exchange interaction 
and the so-called correlated-hopping terms as they 
are of much smaller amplitude at these distances. Note 
also that the Mott criterion at the transition [Ml 123 takes 
in the present situation the form asn^c’^ = {aa)~^ « 0.45 
at the transition, not too far from the value 0.25 for 3- 
dimensional (d = 3) systems. 

4-. Outlook We have described in a rigorous manner the 
behavior of electrons in the short H 2 chains connected 
via the Coulomb interaction with the further molecules, 
as well as utilizing the periodic boundary conditions, 
emulating together an extended-chain behavior. Amaz¬ 


ingly, the results obtained meet some of the features 
one can expect for a three-dimensional hydrogen molec¬ 
ular quasiatomic transformation. Obviously, no de¬ 
tailed realistic phase diagram can be obtained, as both 
the structural details and the effect of zero-point motion 
are still missing. However, a unique connection to the 
Mott-Hubbard transition is established. A direct proof of 
metallicity of the quasiatomic state would involve the de¬ 
termination of optical conductivity |28| . The model for¬ 
mulated here may form a basis for quantitative approach 
to solid hydrogen with a proper account of electronic cor¬ 
relations. In such approach a three-dimensional realistic 
model must be constructed. We think that tackling of 
this problem is possible on a small scale within our Exact 
Diagonalization Ab Initio (EDABI) approach employed 
here and to other problems |29| involving electronic cor¬ 
relations as an essential aspect of their physical proper¬ 
ties. For 3d system, approximate second-quantization di¬ 
agonalization schemes such as the statistically-consistent 
Gutzwiller approximation may be used with concomitant 
single-particle wave-function optimization m- These 
methods can be tested in the present exact limit. 

One should also note that the principal assumption 
made at this point is that the protons form a lattice 
even in the atomic phase, not a proton-electron quantum- 
liquid plasma nniEOHsi]. Such problem is possible to 
tackle within the present model by assessing the optimal 
ground-state energy for a random choice of the proton 
positions. But first, an accurate estimate of the zero- 
point motion amplitude and the corresponding energy 
M must be carried out for an extended system with 
correlations |33j . 

Computer animation of the molecular quasiatomic 
transition is attached as Supplemental material |34| . 

5. Methodological remark The correlated ladder systems 
modeled by the (extended) Hubbard model have been 
studied extensively in the past (cf. e.g. |3S])- Those 
studies differ from the problem discussed here. We start 
with the physical H 2 spin-singlet-insulator chain. Under 
pressure, the spin-singlet states transform into a quasi¬ 
atomic correlated state. Most importantly, we treat the 
correlations within the cluster of n molecules rigorously 
and at each step determine the ladder rung size (bond 
length R) and the equilibrium lattice constant a. Such 
analysis has not been carried out so far and provides also 
the values of the microscopic parameters. 
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